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Abstract 



We study the isospin-singlet neutron-proton pairing in bulk nuclear matter 
as a function of density and isospin asymmetry within the BCS formalism. 
In the high-density, weak-coupling regime the neutron-proton paired state is 
strongly suppressed by a minor neutron excess. As the system is diluted, the 
BCS state with large, overlapping Cooper pairs evolves smoothly into a Bose- 
Einstein condensate of tightly bound neutron-proton pairs (deuterons). In the 
resulting low-density system a neutron excess is ineffective in quenching the 
pair correlations because of the large spatial separation of the deuterons and 
neutrons. As a result, the Bose-Einstein condensation of deuterons is weakly 
affected by an additional gas of free neutrons even at very large asymmetries. 

PACS: 21.65.+f, 74.20.Fg, 24.10.Cn 
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I. INTRODUCTION 



The crossover from BCS superconductivity to Bose-Einstein condensation (BEC) man- 
ifests itself in fermionic systems with attractive interactions whenever either the density is 
decreased and/or the interaction strength in the system is increased sufficiently The tran- 
sition from large overlapping Cooper pairs to tightly bound non-overlapping bosons can be 
described entirely within the ordinary BCS theory, if the effects of fluctuations are ignored 
(mean- field approximation). Indeed in the free-space limit the gap equation reduces to the 
Schrodinger equation for bound pairs. This type of transition has been studied initially in 
the context of ordinary superconductors excitonic superconductivity in semiconductors 
0, and, at finite temperature, in an attractive fermion gas ||. Although the BCS and BEC 
limits are physically quite different, the transition between them was found smooth within 
the ordinary BCS theory. 

More recently it was argued that a similar situation should occur in symmetric 

nuclear matter, where neutron-proton (np) pairing undergoes a smooth transition from a 
state of np Cooper pairs at higher densities to a gas of Bose-condensed deuterons, when the 
nucleon density is reduced to extremely low values. At the same time the chemical potential 
evolves from positive values to negative ones (disregarding a mean field), approaching half 
of the deuteron binding energy in the zero-density limit. This transition may be relevant, 
and could give valuable information on np correlations, in low-density nuclear systems like 
the surface of nuclei, expanding nuclear matter from heavy ion collisions, collapsing stars, 
etc. 

The np pairing effect is largest in the isospin symmetric systems. However, most of the 
systems of interest are, to some extent, isospin asymmetric, i.e., there is usually a certain 
excess of neutrons over protons. In this case the pairing is suppressed, since for neutrons 
and protons lying on different Fermi surfaces the phase space overlap decreases as these 
surfaces are pushed further apart by isospin asymmetry f|,|7|||. Note that this situation 
is quite analogous to pairing of, e.g., spin-polarized electrons. As is well known, an extra 
amount of spin-up electrons over the spin-down electrons is very efficient in destroying the 
superfluidity in the weak-coupling regime. 

This situation is well known also from odd nuclei or nuclei with quasi-particle excitations, 
where the so-called blocked gap equation emerges ||. Indeed, the ground state wave function 
for a np BCS state with N excess neutrons can be written in the following way ||, 



where and are the quasiparticle and bare particle creation operators, and |BCS) is 
the superfluid ground state in the symmetric case. From Eq. ([[]) we see that the presence 
of the extra neutrons entails a suppression of the Cooper pairs in the np BCS state, which 
have the same quantum numbers as the extra neutrons. The variational minimization of 
the Hamiltonian then leads to the so-called blocked gap equation ||, where the window of 
ki, . . . , fc^v states is missing in the integral equation for the gap. Since this window is situated 
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close to the Fermi energy, where most of the pairing correlations are built, the suppression 
mechanism is extremely efficient. This is the case for positive chemical potential. 

The situation is, however, more complex when the chemical potential approaches zero 
and becomes negative, i.e., in the limit of strong coupling. Clearly the np pairs will evolve 
to tightly bound deuterons in this limit, and, because of the isospin asymmetry, they will 
coexist with a dilute gas of excess neutrons. It is physically conceivable that in this limit, 
where the deuterons as well as the extra neutrons are extremely dilute, the Pauli blocking 
effect of neutrons on the deuterons will become negligeable. This may have interesting 
consequences, for example, in the far tail of nuclei, where a deuteron condensate may exist 
in spite of the fact that there the density can be quite asymmetric. Apart from finite nuclei, 
similar physical effects could play an important role in other low-density asymmetric nuclear 
systems. 

The aim of this work is, therefore, a detailed study of the behavior of np pairing in 
asymmetric nuclear matter as a function of density and asymmetry. The paper is organized 
as follows. In section we set up the general frame in terms of the Gorkov formalism at 
finite temperature. In section [TT1| we specify the relevant equations to zero temperature and 
discuss the analytical limiting case in the weak-coupling approximation. Numerical results 
are presented and discussed in section [TV]. Our conclusions are summarized in section M. 



II. BASIC EQUATIONS 

In this section we briefly review the treatment of an isospin-singlet superfluid Fermi sys- 
tem within the Green's function formalism, see also Refs. |5|-§] . We use the standard method 
of Green's functions at finite temperatures, see, e.g., Refs. fl0|-|l4"|| . In this formalism inter- 



acting superfluid systems are described in terms of the Gorkov equations, which generalize 
the Dyson equation for normal Fermi systems by doubling the number of propagators. For 
a homogeneous system the normal and anomalous propagators, G and F, are defined as 

G Ta , T , a ,(k,t-t') = ^(T[iP Ta {k,t)iPl ltJ ,{k,l/)}) , (2a) 

F Ta y a/ (k,t-t') = t(T[Vw(MVv<x'(M')]> , (2b) 

where a — f, J, and r = n,p denote the spin and isospin quantum numbers, respectively. On 
introducing the Matsubara frequencies to v = (2z/+ l)7rT, where T = 1/(3 is the temperature, 
v — 0, 1, 2, . . ., the propagators are conveniently written in the Fourier representation as 

G ra y a/ (k, t) = ^Yl e ~^" tG ra,r'a' (fc, 0J V ) , (3a) 
v 

F T0VT v (k, *) = 4 £ e ~ W " tF ™y«' <»v) ■ (3b) 

V 

Using the short-hand notation D, A, G, and F, for 4x4 matrices in spin-isospin space, the 
Gorkov equations can be written as 



D(fc,u;„) A(fc) \(G(k,u v )\ (l 
At(fe) -D(fc, -u„) J 1 Ft(fe, u v )J 1 



(4) 
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where 



(iu v -£ n ^(k) \ 

iuo u — E n i(k) 

ioj u — £ P f(fe) 

\ iuj v - e pl (k) J 



(5) 



and 1 denotes the four-dimensional unit matrix. Since we are mostly interested in the 
pairing properties at very low density, where the nuclear mean field plays a minor role, the 
single-particle spectrum adopted in this work is simply the kinetic one, 



k 2 

e n ^{k) = e ni (k) = e n (k) = /i n = e k + 5/j 

2m 

k 2 

£ pT (fc) = e pl (k) = e p (k) = — - n P = e k - 5n 



(6a) 
(6b) 



where e k = k 2 /2m — fx, pL — (fi n + yU p )/2 is the average chemical potential between protons 
and neutrons, and 5fi = (fi n — fi p )/2 is the associated shift. From now on we assume nuclear 
matter with neutron excess so that 5fi > 0. 

In this article we concentrate on pairing in the dominant isospin-singlet channel that can 
be described by a gap function A(fc) with the structure 



/ A + iAi \ 

-A + zAi 

-A - zAi 
\A -iAi 



(7) 



which is a particular case of a unitary state P,|TR|IB 



A f A = A 2 l, A 2 = v / detA = A 2 + A? 



(8) 



It allows in principle the coexistence of spin-singlet (S = 0) and spin-triplet (S = 1) pairing 
correlations Ag. In the low-density region that we are interested in, pairing is realized in 
the spin-triplet s-wave channel 3 SD\, i.e., A = iA±. 

The anomalous propagator has the same spin-isospin structure as A, whereas the 
normal propagator G is diagonal in the spin-isospin indices. Taking this into account, the 
system Eq. (||) can be inverted with the solution 



G n (k,u u ) 
G p {k,u v ) 



iu v + £ p (fe) 



(iu v - E+) {iu v + E k ) 
iu v + e n (k) 



E k ) + E+) ' 

-A(fc) 
EH) K + E+) ' 

+A(fe) 



{iuj v - E+) {iu u + E k ) 



(9a) 
(9b) 
(9c) 
(9d) 
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where 



E£ = E k ± Sfi = J e 2 k + A(k) 2 ± Sfi . 



(10) 



The isospin asymmetry thus lifts the degeneracy of the quasiparticle spectra, leading to two 
separate branches for protons and neutrons. This is sketched in Fig. [I]. In the region around 
/i, where E k < <5/i, the superffuid state is unstable against the normal state and indeed it 
accomodates the unpaired neutrons. This "isospin window" is responsible for the pairing 
gap suppression in the density domain where nuclear matter is superfluid. It corresponds to 
the blocking effect for the pairing in nuclei [[J. 

After analytical continuation of the Green's functions in the complex cj-plane one calcu- 
lates the spectral function A(k, u), which is given by the discontinuity of G across the real 
axis, 



A T (k,u) = 2tt 



(ii) 



and then the density of particles 

duo 



Pr 



2 E 

k 

2 £ 



(2*) 



A T (k,u)f(u) 



.1 + — 
2 V E k 



E k 



(12) 



where Ylk = I d 3 k/(2ir) 3 and f(E) = [1 + exp(/3E)] 1 is the Fermi distribution function. 
Analogously, the discontinuity of the anomalous propagator F\ 



B(k, uo) = 2tt [S(u - E+) - 5{uo + E t 
yields the gap equation 

A(k')=J2(k'\V\k) f^B(k, U )f{. 



-M A(fe) 



2E k 



(13) 



Affe 



(14) 



Using the angle- averaging procedure, which is an adequate approximation for the present 
purpose (see Ref. |I6[), the BCS gap equation for asymmetric nuclear matter can be derived: 



(15) 



M*0 = - EE 7 «'(*' fc ')^ I 1 - ^ W - TO! ' (M' = o,2) 



where E\ = el + A(k) 2 and A(k) 2 = A Q (k) 2 + A 2 {k) 2 is the angle-averaged neutron-proton 
gap function. The driving term, Vu>, is the bare interaction in the 3 SDi channel. We use 
for the numerical computations in this work the Argonne Vu potential fTi 
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From Eq. (|T2|) the total density p = p n + p p and the neutron excess 8p = p n — p p can 
easily be derived and one gets 



p = 2^Tn fc , n k = 1 - Jj- [1 - /(£+) - /(£,")] , (16a) 
k k 

6p = 2j2^ k , 6n k = f(E-)-f(E+). (16b) 



For a convenient parametrization of the total density, we also introduce the quantity kp = 
(37r 2 p/2) 1//3 , which is however, apart from the isospin symmetric system, not to be identified 
with a Fermi momentum. 

Introducing the anomalous density 



Mk) = « fc <-^ = ^ [1 " /(^ + ) " > ( 17 ) 



and making use of Eq. ( |16a| ), the gap equation, Eq. (|15l), can be recast in the Schrodinger-like 
form 

—ipl{k) + (1 - n k ) J2 J2 V "'^ k ')Mk') = 2/i Mk) ■ (18) 
y v 

In the limit of vanishing density, n k — > 0, this equation goes over into the Schrodinger 
equation for the deuteron bound state 0. The chemical potential 2p = p p + p n plays then 
the role of the energy eigenvalue. 

Let us finally remark that in this article we do not study the competition between and 
possible coexistence of isospin singlet and triplet pairing |Tj|. It is clear that T — 1 nn 



and pp pairing gaps, which are unaffected by the isospin asymmetry, will be larger than the 
T = gap at sufficiently large asymmetry. Our results are thus valid up to the asymmetries 
(yet unspecified) at which the isospin triplet pairing becomes dominant. On the other hand, 
as we show below, for very low density systems the suppression mechanism due to the shift of 
Fermi surfaces is ineffective and it is reasonable to assume that the np pairing will dominate 
other channels in a wide range of values of the asymmetry. 

We also restrict ourselves to work with a free single-particle spectrum, which is adequate 
at low density. At higher density, when the effective mass deviates substantially from the 
bare mass, renormalization of the single-particle spectrum should be taken into account 
This affects the magnitude of the pairing gap in general and the critical asym- 



metries at which the pairing effect disappears ||. In view of this approximation, the results 
shown below should be considered as only qualitative in the high-density region. 



III. ZERO TEMPERATURE GAP EQUATION 

In order to disentangle the isospin effects from the thermal ones, we focus in the following 
on the limit of zero temperature, where f(E£) = and f(E^) = 1 — 9(E^), 6 being the 
step function. In this limit Eqs. (|15D and (0) reduce to the following ones, 
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M*0 = " E E V ^ k ')^ 6 ( E k) . (19a) 



P = 2 E| 1 -# fl K) > (19b) 



5p = 2E[l-^)] • (19c) 

In general these three coupled nonlinear equations have to be solved numerically, maintain- 
ing the self-consistency. Before presenting the numerical results, let us however first discuss 
in some detail the physical content of these equations. The physical interpretation of the 
formalism is that in asymmetric nuclear matter a superfluid state of np Cooper pairs with 
density p — 5p = 2p p coexists with a gas of free neutrons with density 5 p. Since this latter 
is occupying a region around the Fermi surface [determined by 1 — 6(E^)}, the momentum 
space available for the pairing is reduced. Thus the effect of isospin asymmetry is to reduce 
the magnitude of the energy gap; with increasing neutron excess the superfluidity rapidly 
disappears @,0,§. 

The solution of the gap equation, Eq. (0), assumes different properties according to the 
value of the chemical potential. One may distinguish three domains: (i) /j > A (weak- 
coupling, pairing regime), (ii) p ~ (Mott transition), and (iii) p < (strong-coupling, 
bound-state regime). Approximate analytical results can be found in the cases (i) and (iii), 
and are described in the following: 



(i) Weak-coupling regime 

In the region p ^> A the gap equation in the form of Eq. ( |19a|) can be solved in the usual 
weak-coupling approximation [21], in order to gain some insight in the qualitative behavior 



of the pairing properties in the asymmetric case. This approximation is not very accurate, 
but it preserves the physical content of the exact solution. 

First, one notes from Eq. (|19cj ) that the unpaired neutrons are concentrated in the energy 
interval [p — Se, p + Se], with the half-width 



Se = ^Sp 2 - A 2 . (20) 

This interval is free of protons and does not contribute to the pairing interaction, whereas 
outside of it neutron and proton distributions are equal and given by the BCS result, see 
Eq. ( |19b| ). This situation is sketched in Fig. ^. One notes at this point that the condition 
5p > A has to be fulfilled in order to generate an asymmetry. Indeed for A, Se <C p the 
momentum distributions of the two species are very sharp (practically Fermi distributions) 
and one obtains as relation between the asymmetry a and Se: 

a p n - Pp ^ {fo +Se -fo' 6£ ) d ^ ^3Se 



7 



i.e., the asymmetry is directly proportional to 5s. Thus we see that 5s has the interpretation 
of an "effective" difference of chemical potentials, since it is 5s and not 5fx which determines 
the density of excess neutrons for A ^ 0. Only in the normal system one has 5s = 5li, 
whereas pairing screens 8fi such that 8e < 5fx, and it needs a critical finite 8fi = A in order 
to obtain a finite density of excess neutrons. 

Coming now to the gap equation, with a contact interaction V and energy cutoff s c , 
Eq. ( |19a| ) reads 

{2mf/ 2 V [** I s + li 



^ 2 Vl ' ' V s 2 + A 2 



tem-Sehhrrh- ( 22 ) 



For A C /i < e c , the integral can be approximately calculated, yielding 

8s + 8fi 



N(0)Vln 



2JJx~S c 



(23) 



where N(0) = mkp^-n 2 is the level density at the Fermi surface of one type of particles. 
Therefore 

6fi + 6e = const. = A ^> A 2 = A 2 - 2A 5e , (24) 

where Ao is the value of the gap in the symmetric system of the same density. We again see 
that it is 8e, and not 5ii, which determines the reduction of A from its symmetric value. The 
gap decreases very rapidly and disappears when the width of the window reaches 25 e = A , 
i.e., the size of the gap at symmetry. This is shown in Fig. H(a). 



We can now combine Eqs. (|2l"D and (24) in order to specify the dependence of the gap 
on the asymmetry: 



The superfluidity vanishes smoothly, but with an infinite slope at 

"max = -J— , (26) 

which in the weak-coupling limit is a very small value. 

We can also determine the dependence of the gap on the difference of chemical potentials 
8 pL [0,g2|: 



A ^ 1. (27) 



A V A 

In the symmetric system one has 5fi = A , 8e = 0, while increasing the asymmetry the gap 
decreases with decreasing 8fi and increasing 8s and vanishes at 8fi = 8s — A /2. 

Physically, this behavior can be understood by noting that the chemical potential dif- 
ference 25 fx = fx n — Hp is the energy one must invest in order to remove a proton from the 
system and then to insert a neutron instead, in other words the symmetry energy. However, 
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in order to do so in the superfruid system with a neutron excess (even very small), one must 
necessarily break a pair in order to remove the proton, but then one does not gain energy 
by adding the neutron to the top of its Fermi sea. That is why Sfi > A in the pairing regime 
and 5p, > \J pi 2 + A 2 in the bound state regime (the equalities holding in the symmetric 
system). So, clearly in a superfruid system of any density and asymmetry, 5/j, can never be 
zero, because it involves breaking a pair (or bound state). 

Then, successively increasing the asymmetry by breaking pairs in this way and filling up 
the Fermi sea of excess neutrons, the pairing correlations are destroyed, and consequently one 
must invest less energy in order to exchange p <-> n. Therefore Spt decreases with increasing 
asymmetry (while Se increases), up to the point where, at 5fi = 5e = A /2, the np gap 
disappears completely and 25fi is simply given by the difference of Fermi energies of the 
two noninteracting Fermi gases at this asymmetry. The system is then in the normal phase, 
5e = 5fi, which means that in order to further increase the asymmetry the parameter Sp, 
has to increase again, up to 8u = u, corresponding to pure neutron matter. Consequently, 
values of Sp, between A /2 and A are present in the superfluid as well as in the normal 
phase, corresponding to different compositions of small and large asymmetry, respectively. 
The relations between A, 5fi, and 5e are sketched in Fig. |. 

We comment finally on the analogy to an electron system in a magnetic field. Whereas 
in the case of nuclear matter discussed before, an increasing asymmetry is imposed on the 
system, leading to a smooth disappearance of the gap at a certain maximum asymmetry, the 
situation for the electron system is qualitatively different: Varying in this case the magnetic 
field is equivalent to imposing the chemical potential difference 8p, instead of the asymmetry, 
which now corresponds to the magnetization of the system. Therefore one observes in this 
case a first-order phase transition at 5u = A , when the pairs are suddenly broken up and the 
system jumps immediately from the symmetric superfluid phase to the magnetized normal 
phase. A magnetized superfluid phase is thermodynamically unstable in this case. 

(ii) Mott-transition regime 

Let us now discuss the situation when the density p decreases. As we have shown in 
H, simultaneously the chemical potential p, decreases (disregarding the Hartree-Fock shift) 
and at a certain (very low) density (of the order of p ~ po/100, where po ~ 0.17 fm~ 3 is 
the nuclear matter saturation density) p passes through zero, corresponding to the Mott 
transition of the deuteron |H|. This also remains true in the asymmetric case; in fact, the 
weak-coupling approximation used above predicts that in pure neutron matter this happens 
at ji n = (7r/2) 2 /i£>, where ud = Ed/2 ~ 1.1 MeV is half of the deuteron binding energy. 
This corresponds to a density p ~ 0.0016 fm~ 3 . 

It seems clear that the Pauli blocking starts to loose its efficiency once the left-hand 
side of the window [p — 5e, p + 8s\ passes beyond zero, because then only part of the 
window participates in the blocking. When p has reached negative values, less than half of 
the blocking window is actually effective. While in the symmetric case nothing particular 
happens when the np Cooper pairs change their character to bound deuterons in the low- 
density limit, the asymmetric case needs special attention. Here both <5p and p vary strongly 
with asymmetry (at fixed total density), and the results have to be found numerically, solving 
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the coupled Eqs. ([19|). 



(Hi) Strong-coupling regime 

Once the density p becomes so low that simultaneously the mean distance among the 
deuterons as well as among deuterons and excess neutrons becomes much larger than the 
deuteron radius, the excess neutrons can not excert any significant influence on the deuteron 
wave function. This is quite opposite to the weak-coupling case p ^> A considered above, 
where we have seen that only a slight neutron excess destroys the np Cooper pairs. 

For negative chemical potentials Eq. (^) is not valid, since the Fermi surface drops 
into the unphysical region. However, the low-density limit p — > of the BCS equations 
provides simple analytical expressions, once the density is so low that the chemical potential 
is close to the asymptotic value p — > —po = —Ed/2 m —1.1 MeV. In some sense now the 
gap equation (|19a|) and the density equations ( |19b| , |19c|) interchange their roles, because the 



gap equation goes over into the Schrodinger equation, determining the chemical potential, 
whereas the value of the gap can be extracted from the equations for the density, in the 
following manner: The density of protons is given by 



(2m) 3 / 2 f°° , ^ 1 
Pp = 0Z2 / de Ve 



2 71 " -'e min <|/i| 2 



e - p 



V(e-/x)2 + A2. 

The numerical value of the last integral being tt/2, one obtains for the asymmetry 



with 



a = = 1 - ^ (29) 

P P Pd 



(2mp 



_ K „.„rDj — _ 0.00049 fm" 3 . (30) 

07T 



1 3/2 



The low-density result is therefore |23 



A (fta) = Jee-°>, pi) 

Pd y P m 

and thus a gap exists for any asymmetry. 

The ineffectiveness of the neutrons in the p — > limit also becomes evident when consid- 
ering Eq. (0). Indeed in the limit p n , p p —>■ also the momentum distribution vanishes 
and then the equation goes over into the free deuteron Schrodinger equation, independent 
of the asymmetry of the system. Since the Pauli blocking becomes less and less important 
as the density decreases, the asymmetry can become larger without destroying the super- 
fluidity. As we mentioned already in the introduction, this aspect can become important in 
the far tail of the nuclear densities or in other low-density nuclear systems. 
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In the presence of a neutron excess the Bose-condensed deuterons occupy the negative 
energy state E pa —Ed and the free neutrons the positive energy states. With increasing 
asymmetry the neutrons accomodate in the next positive energy states according to the 
Pauli principle. As we see from the wave function in the introduction, even in the low- 
density limit the excess neutrons stay antisymmetrized with the neutrons bound in the 
deuterons. Therefore one cannot distinguish between bound and unbound neutrons, the 
chemical potential of the neutrons is always the one of the unbound ones which is tending 
to zero as p approaches zero. On the other hand the proton chemical potential p p tends 
to the binding energy of the deuteron, since it is the binding energy of the system per half 
the number of particles bound in the deuterons. Therefore the mean chemical potential 
P — + A t p)/2 tends to half the binding energy of the deuteron such that the eigenvalue 
2p of Eq. (|IBD hits precisely the eigenvalue of the deuteron at p = 0. 

To summarize, the behavior of the various chemical potentials in the low-density limit 
is p n — > 0, p p — > — E D , p — > —Ed/2, 5p — > E D /2, independent of the asymmetry! 



IV. NUMERICAL RESULTS 

We discuss now the results that were obtained by numerical solution of the system of 
equations (|19D, using the Argonne V u potential as the bare interaction. 

We begin, in Fig. f|, with an overview of the resulting gap in the density-asymmetry plane. 
As discussed in section [TTT], in the high-density pairing regime the superfluidity vanishes at 
a finite asymmetry a max (p), decreasing with increasing density, whereas in the low-density 
bound-state region (below kp pa 0.35 fm -1 , p ~ po/100) a gap exists for any asymmetry, 
roughly following the analytical result, Eq. (|31|). 

In Fig. [|, we show the momentum distributions, gap functions, and anomalous densities 
at three different densities p = 10~ 4 , 10~ 2 , 10 _1 fm~ 3 , representative of the low, intermediate, 



and high-density regions mentioned in section |TJ. The top panels of the figure show the 
momentum distributions of the superfluid protons. As discussed before, see Fig. |2], there is 
a proton-free region that is at high density centered around the Fermi momentum, and at 
low density excludes all low-momentum states, representing the Fermi sea of non-superfluid 
excess neutrons. The resulting variation of the gap function in the relevant regions of 
momentum space is rather small at low density, however, it is competing with the kinetic 
energy in the self-consistent determination of the normal region, which is determined by 
the variation of E k = y ' e\ + A 2 . The anomalous density develops a characteristic peak at 
the Fermi momentum only in the high-density regime, whereas at low density it exhibits a 
smooth variation typical of the deuteron wave function. 

In Fig. |^, we show in more detail the gap in symmetric matter, together with the an- 
alytical approximation Eq. (|3l|) (top panel); the maximum value of asymmetry at which 
the superfluidity disappears, in comparison with the analytical estimate Eq. (|26|) (mid- 
dle panel); as well as, in the bottom panel, the variation of the various chemical poten- 
tials along the line [p, a max (p)]. The asymptotic behavior predicted in the previous section, 
(6p, p n , p, p p ) — > (1, 0, —1, —2)p D , is observed for p — > 0. 

We finally discuss briefly the situation at finite temperature, that was treated in detail 
in Refs. H&fl. As long as T < A, the temperature effects are small. However, in the 
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low-density limit, when A — > 0, a qualitatively different behavior from the zero temperature 
case can be observed: In Fig. ||| we have seen that in this limit /x„ — and p p — > — 2/zq, 
due to the coexistence of the deuterons with a Fermi sea of free neutrons, see the discussion 
in section |l|(i). At finite temperature, however, part of the deuterons are broken up and 
therefore this argument does not apply any more. This can be seen from Eq. ( |16b| ) for Sp. In 
the limit p, 5p —>■ we have A = 0, since at finite temperature the system becomes normal 
at a certain critical density. Therefore Eq. ( |16b| ) reads 

= 2 ]T [f(e k - 5p) - f(e k + 5p)} (32) 

k 

and one clearly sees that this relation can only be fulfilled, at any finite temperature, for 
5p = 0, i.e., p n = p p . On the other hand, we see from Eq. ( |18D that p — > —po- Therefore 
we have for T / the limit p,p n ,p p — ► —po- This behavior is clearly born out from 
the numerical calculation, Fig. |7], which shows the density dependence of p n and p p for 
fixed asymmetry a = 0.1 and various temperatures T = 0.1,0.5,1.0 MeV. It is seen how 
with vanishing temperature the situation shown in the lower panel of Fig. ^ (for a = 1) is 
approached. 

V. CONCLUSIONS 

In this article we extended a previous study [|J of the transition from a neutron-proton 
BCS superconducting state to a Bose condensate of deuterons in symmetric nuclear matter 
to systems with isospin asymmetry. This is an important aspect, since most of the low- 
density nuclear systems, such as tails of nuclear density distributions in nuclei, have a strong 
neutron excess. In the high-density weak-coupling regime, p ^> A, even a small asymmetry 
is sufficient to suppress the pairing correlations completely via the Pauli blocking effect. 
However, we find that in the situation where the density drops so low that the chemical 
potential turns negative and deuterons start to form, the isospin asymmetry is much less 
important, i.e., the system supports pair correlations (in the form of a deuteron condensate) 
for much larger asymmetries. 

Indeed one can argue that, once the density is so low that the spatial separation between 
deuterons and between deuterons and extra neutrons is large, the Pauli principle is ineffective 
and, hence, asymmetry can not destroy any longer the binding of neutron-proton pairs. 
Our numerical calculations fully confirm this conjecture. Most easily this effect can be 
understood by examining the phase space blocking effect on the anomalous density. In the 
high-density regime the latter quantity shows a peak structure as a function of the mean 
chemical potential of neutrons and protons p. Once the asymmetry is imposed the Pauli 
blocking effectively cuts out the major part of this peak structure. When the chemical 
potential turns negative, there is no peak any more in the relevant physical region of the 
anomalous density and the Pauli blocking looses its efficiency, which enables the proton- 
neutron pairs to condense in the very low density regime even in the presence of a large 
neutron excess. We anticipate that this aspect may be important for the understanding of 
the far tails of density profiles of exotic nuclei, the expanding (asymmetric) nuclear matter 
in heavy ion collisions, and other low-density nuclear systems. 
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FIGURES 
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FIG. 1. Sketch of the quasiparticle spectrum. For e < // the lower (upper) solid line correponds 
to neutrons (protons) in the normal state. For e > fi the lower (upper) solid line corresponds to 
protons (neutrons) in the normal state. The corresponding dashed lines refer to nucleons in the 
superfluid state. 
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FIG. 2. Momentum distributions of neutrons and protons in asymmetric superfluid matter. 
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FIG. 3. Dependence of gap A (a) and chemical potential difference 5/j, (b) on the width of the 
blocking interval Se within the weak-coupling approximation, see Eq. (p^). The solid lines denote 
the superfluid phase and the dashed lines the normal one. 
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FIG. 4. The pairing gap as a function of total density [kp = (3ir 2 p/2) 1 / 3 ] and asymmetry. 
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FIG. 5. Proton momentum distributions n p (k), gap functions A(k), and anomalous densities 
ip(k), for different total densities p = 10~ 4 , 1CP 2 , lO^ 1 fm~ 3 , and different asymmetries. The dots 
indicate the boundaries of the interval containing the neutron excess in each case. 
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FIG. 6. Top panel: The 3 SDi gap in symmetric nuclear matter as a function of total density p 
or equivalent Fermi momentum = (3ir 2 p/2) 1 / 3 (solid line). The dashed line shows the analytical 
approximation Eq. ( pip . Central panel: The maximum asymmetry at which a gap exists (solid line) 
and the approximation Eq. ( |26| ) (dashed line). Lower panel: The values of the various chemical 
potentials p, Sp, fj, n , p p , corresponding to the asymmetry a max displayed in the panel above. 
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FIG. 7. Chemical potentials of neutrons (upper curves) and protons (lower curves) at a = 0.1 
as functions of density for different temperatures. 
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